Preprint typeset using I^TgX style emulateapj v. 5/2/11 
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ABSTRACT 

We investigate how strong lensing of dusty, star-forming galaxies by foreground galaxies can be used 
as a probe of dark matter halo substructure. We find that spatially resolved spectroscopy of lensed 
sources allows dramatic improvements to measurements of lens parameters. In particular we find that 
modeling of the full, three-dimensional (angular position and radial velocity) data can significantly 
facilitate substructure detection, increasing the sensitivity of observables to lower mass subhalos. We 
carry out simulations of lensed dusty sources observed by early ALMA (Cycle 1) and use a Fisher 
matrix analysis to study the parameter degeneracies and mass detection limits of this method. We 
find that, even with conservative assumptions, it is possible to detect galactic dark matter subhalos 
of ~ 10^ Mq with high significance in most lensed DSFGs. Specifically, we find that in typical DSFG 
lenses, there is a ~ 55% probability of detecting a substructure with M > 10^ Mq with more than 
5(7 detection significance in each lens^ if the abundance of substructure is consistent with previous 
lensing results. The full ALMA array, with its significantly enhanced sensitivity and resolution, should 
improve these estimates considerably. Given the sample of ~ 100 lenses provided by surveys like the 
South Pole Telescope, our understanding of dark matter substructure in typical galaxy halos is poised 
to improve dramatically over the next few years. 

Subject headings: dark matter — gravitational lensing — galaxies: dwarf — galaxies: luminosity 
function, mass function — galaxies: structure — 



1. INTRODUCTION 

The infiationary ACDM model has proven spectac- 
ularly successful in describing the observed large-scale 
structure of the universe. Precise measurements of the 



mean expansion history of the universe (e.g. Suzuki et al 
2012), t he anisotropics of the cosmic microwave back- 
ground (Komatsu et al. 2011), and the clustering of 



galaxies at low r edshitt z 5. 1 (e.g. |Reid et al.||2012 
Sanchez et al.| |2012) are all consistent with a simple see 
nario involving a nearly scale-invariant spectrum of cur- 
vature perturbations generated during infiation, whose 
growth over cosmic time is governed by gravitational in- 
stability. 

While measurements of large-scale (> 10 Mpc) struc- 
ture have yielded the most stringent constraints on this 
cosmological model, there is considerable information to 
be gleaned from observations of structure on smaller 
scales as well. For example, the detailed shape of the in- 
fiationary potential in fiuences the sha pe of the primordial 
power spectrum (e.g. Dodelson] |2003 ) , motivating efforts 
to measure the scalar spectral index Uq and its running 
dus/dlogk. For example, a sharp feature in the infia- 
tionary potential will generate features in the primor- 
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dial power spectru m, breaking its near scale-invariance 
( Kamionkowski & Liddle 2000). In addition, the particle 
physics ot dark matter affects the shape of the matter 
transfer functions, leading to a Jeans-like suppression of 
structure below the free-streaming scale of the dark mat- 
ter particles. In WIMP-like scenarios, this free-st reaming 
scale is as small as 1 comoving pc ( jLoeb fc Zaldarriaga] 



2005 ), however in alternative dark m atter models this 
damping scale can be much larger (e.g. [Bode et al. 2001 



'Abaza.iian'2006UCembranos et al.|2QQ5lfKaplinghatl2005, 

Cyr-Racine & Sigurdson 2012^. Accordingly, measure- 
ments of the small-scale power spectrum can in principle 
provide a wealth of information about the physics of the 
early universe. 

To date, the most constraining bounds on the small- 
scale power spectrum have b een derived from observa- 
tions of the Lyman-a forest (Seljak et al. 2006). How- 
ever, it will be difficult for future LyAI^ observations 
to improve significantly on existing bounds because cur- 
rent, high-resolution spectra already have resolution ap- 
proaching the Jeans scale of the intergalactic medium 
at z ~ 2 — 3. Another probe of small-scale structure is 
the abundance of low-mass dark matter halos and subha- 
los. High-resolution N-body simulations 
Stadel et al 



2007 



Diema nd et al 
Navarro et al. ,2010l) have 

rDM. 



I ^ ,2008; Sta del et al.,2009, —\^rr.: 

revealed that dark matter halos in hierarchica!! CJDiVi coS' 
mologies are expected to contain copious substructure, 
primarily in the form of gravitationally self-bound, dy- 
namically cold subhalos with a spectrum of masses. The 
abundance of these subhalos is expected to depend on 
the amplit ude and shape of the sm all-scale power spec- 
trum (e.g |Zentner fc Bullock 2003), meaning that mea- 
surements orTiaIo"^u5striictiire^^ constrain early- 
universe physics. 
Numerous groups have attempted to constrain the 
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abundance of low-mass subhalos by detecting faint, DM- 
dom inated dwarf sa tellite galaxies in the Local Group 
for a review). 



(see Kravtsov 
forts 



2010 



Despite heroic ef- 

the nuniber of detected dwarf satellites falls well 
below the expected abundance of low-mass subhalos for 
a typical galactic halo. This discrepancy has been termed 
the "Satellite Problem" for CDM models, and has led to 
considerable speculation that the physics of dark mat- 
ter might not be described well by the idealized, non- 
interacting, absolutely cold CDM model. However, as- 
trophysical processes can plausibly suppress the star- 
formation efficiency of low-mass halos and subhalos, ren- 
dering many of them invisible to optical surveys. To 
distinguish between astrophysical solutions, and particle- 
physics solutions to the Satellite Problem, a purely grav- 
itational method for detecting potentially dark subhalos 
is required. One possible means of detecting dark sub- 
structure in galactic halos is to search for its dynamical 



effects on c old stellar streams ( Carlberg et al^2011 Yoon 
et al.||2011 ). Another gravitational probe ot dark matter 



substructure is strong gravitational lensing, which is the 
focus of this work. 



|Mao & Schneider) ( [1998 ) were the first to suggest that 
strong gravitational lenses could be used to detect dark 
matter substructure through anomalous flux ratios of 
multiply imaged quasars. Many groups subsequently fol- 
lowe d up on this idea, b oth theoreticall y and obs ervation- 
ally dDalal fc Kochanek 2002; Kochanek & Dalai 2004; 
Chiba et al. 2005; Minezaki et al. 2009; MacLeod et al. 
2009: Miranda & Maccio 2007: Chen et al. 2007, 2011, 



Rozo et al. 1^ 007; Keeton eTliL]|20061 |Chen et aL 



2003 



Moustakas fc Metcalf 2003; 'Metcalf et al."2004 | [Metcalf 
2002; Schechter fc Wambsganss 2002; Koopmans et al. 



■ 2002, Metcalf fc Madau 2001; Keeton et al. 2003; Fadely 
FKeeton 2011; Keeton et al. 2005; Xu et al. 2009; Ve; 



Keeton :^Uii; Keeton et al. :^UU5; Au et al. :^UU9; Veg- 
etti et al.fm2f .ln particular, Dal ai fc KochaiieE[^2"QQ2f 



analyzed a sample of qu adruply imaged rad io quasars 
from the CLASS survey ( [Browne et al.|[2QQ3 ), and from 
the preponderance of flux anomalies in these lenses, con- 
cluded that a large fraction (~ 1 — 2%) of the projected 
mass at the Einstein radius, ^ 5 kpc, must be in the 
form of local substructure. The uncertainties on this 
measurement were quite large, however, due to the small 
sample size used in that analysis. Subsequent work has 
not significantly improved the bounds on substructure 
from strong lensing, mainly due to the difficulties in using 

ftical lenses to study substructure (Kochanek fc Dalal| 
04), however Vegetti et al. (2012) report a recent sub- 
lo detection using an extended optical galaxy-galaxy 
lens system. 

Fortunately, in recent years a new class of lensed 
sources has been detected: dusty, st ar-forming ga laxies 
(DSFGs) at high redshift, z - 2 - 5 ( |Vieira et al.pOlOt 
[Negrelloet al.|201 0). DSFGs are a class of luminous and 
prodigiously star-forming galaxies located at high red- 
shift (z > 1). They are enshrouded in dust (Blain et al 



2002 



Lagache et al. 2005 ) and conta i n massive reservoir s 



of molecular gas (Greve et al. 2005; Carilli et al^|2011 ). 



The molecular gas in these galaxies is excited by the in- 
tense high energy emission of the active star forming re- 
gions. Most commonly, CO, H2O, HCN, and H CO+ are 
the molecular lines ob served in these galaxies ( Solomon] 
|& Vanden Bout 20051). In addition to molecular lines, 
atomic fine structure lines such as ionized carbon ([Cii]) 



has long been known to be a dominant cooling mecha- 
nism for DSFGs, and in some cases can account for 10~ ^ 
of their total infrared luminosity (iStacey et al. 1991). 
The intrinsic sizes and morphologie s of DSl^'Gs are not 
well understood, but some st udies ( Chapman et al 



Tacconi et al 

suggest typical radii of 



2006 Young er et al. 20 08; Engel et al 
1 kpc tor these galaxies. 




submm emission from DSFGs is believed to be dominated 
by multiple compact clumps of intense star formation 
spread out over t he ex tent of the galaxy. For example, 
"Swinbank et al. ( 2010| ) derive upper limits of ~200 pc 
for the diameters (l^wHM) of star forming clumps in a 
lensed DSFG. 
The bri ghtest DSFGs were predicted to be strongly 



lensed (Blain 1996 Negrello et al. 2007 Hezaveh fc 



Holder 20lT|r Tins prediction was contirrned by high- 
resolution follow-up imaging of the brightest extragalac- 
tic sources found by the South Pole Telescope (SPT) 



iCarls trom et al 



2011 Vieira et al. 2012 



and the hersc nel Space Ubservatory ( [iNegrello 



Hezaveh et al 



et al. 2010). Recently Hezaveh et al. (2012) presented 



lens models for four lensed DSl^'Gs observed with ALMA 
in Cycle 0. Even with very short total observing times 
(~ 20 s) and only ~ 15 antennas, they were able to de- 
rive tight constraints on the mass distributions of the 
foreground lenses, thanks to the extreme brightness of 
the sources combined with ALMA's high sensitivity. The 
stringent constraints that are possible using snapshot ob- 
servations with the small subset of telescopes available in 
Cycle strongly suggests that the full ALMA array could 
provide an unprecedented view of the detailed mass dis- 
tributions of lensing galaxies. In particular, ALMA ob- 
servations of lensed DSFGs could revolutionize the study 
of dark matter halo substructure in lens galaxies. 

Lensed DSFGs are a particularly attractive popula- 
tion for substructure studies, due to their great abun- 
dance (compared to quasar lenses). In addition, since 
almost all of the UV and optical emission in the DS- 
FGs is absorbed and reradiated by dust at longer wave- 
lengths, these sources are almost completely invisible in 
optical images. If substructure is detected gravitation- 
ally, deep optical imaging would place stronger limits 
on the mass to light ratios of the galactic satellites than 
would be possible if the source galaxies were bright in op- 
tical bands. Moreover, the high redshifts of the sources 
permit a wide range of possible lens redshifts, potentially 
allowing constraints on any redshift evolution in the sub- 
structure population. 

In this work we study the feasibility of using ALMA 
to detect dark matter substructure in the halos of lens 
galaxies. In particular, we study how the spatially re- 
solved spectroscopy provided by ALMA allows us to re- 
solve source structures and dramatically increase the sub- 
structure detection sensitivity. In ^we discuss the bene- 
fits of spectroscopically resolved interferometric observa- 
tions of lensed DSFGs. In ^ we describe our simulations 
of ALMA observation of lensed DSFGs. In Q and [5] we 
present results of our simulations, and conclude in q6| 

2. SPATIALLY RESOLVED SPECTROSCOPY OF LENSED 
SOURCES 

One of the most important properties of lensed sys- 
tems which determines the sensitivity of observations to 
substructure lensing is the size of the source. A source 
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that is completely uniform across some scale rgrc will 
be insensitive to lensing perturbations from structures 
on scales small compared to rgrc- When studying dark 
matter substructure, it is therefore advantageous to use 
lensed sources that intrinsically have structure on an- 
gular scales of order milliarcseconds, comparable to the 
Einstein radii of the subhalos of interest. For example, 
radio quasars are excellent sources for studying substruc- 
ture, since their intrinsic sizes {6 < 0.1 mas) are much 
larger than the Einstein radii of individual stars but 
smaller than the Einstein radii of subhalos. In contrast, 
optical QSO's have much smaller angular sizes, mean- 
ing that their lensed images can be affected by stellar 
microlensing which c an be difficult to disentangle from 
substructu re lensing ([Kochanek fc Dalai 20041 Morgan' 
|e^al.|2012| [Sluse et al.|2012). Similarly, typical optically 
bFight galaxies have much larger angular sizes, of order 
~ f , rendering most galaxi es insensitive to substructure 
pert urbations (however s ee Vegetti & Koopmans (2009) 
and Vegetti et al.| ( [2012 ) for a method which uses the 
compact details of the structure in the source for sub- 
structure detection). 

Naively, DSFGs would appear to be ill-suited as 
sources for substructure lensing. Typical DSFGs are be- 
lieved to have physical radii of order 1 kpc, correspond- 
ing to an angular scale of order 0.2^', much larger than 
the Einstein radii of all but the most massive subhalos. 
However, these galaxies do not have a smooth and uni- 
form morphology. They typically are believed to be com- 
posed of multiple compact knots of star formation, based 
both on theoreti cal and observational grounds ( |Swin- 
bank et al. 2010). In principle, these compact clumps 
should be sensitive to lensing perturbations on smaller 
scales than the galaxy as a whole. In practice, however, 
the superposition of a large number of blended, over- 
lapping source clumps becomes indistinguishable from 
an extended source, thereby reducing sensitivity to sub- 
structure lensing. In order to fully exploit the potential 
of lensed DSFGs as a probe of dark matter substructure, 
we require some method to decompose the source emis- 
sion into its constituent compact clumps. 

Fortunately, it is possible to perform such a decom- 
position i n velocity space. The idea is an alogous to the 



metho d of |Moustakas fc Metcalf| ( |2003[ ) and |Metcalf et al 
(|2004|). Spatially resolved spectroscopic observations of 



s show strong velocity gradients, some indicating 
fast -rotating disks, and others indicating major merg- 
ers fEngel et al.||2010[ [Riechers et^|2011[ [Hodge et al. 
'20l2ir If the line-ol-sight velocity ottset between star 



forming clumps is larger than the velocity dispersion of 
single clumps, then observations at different frequencies 
will allow us to distinguish the emission from distinct 
clumps. In other words, the emission in a narrow fre- 
quency window will come from a region in the source 
much smaller than the overall extent of the DSFG. Ef- 
fectively, we can use spectroscopic resolution to enhance 
our spatial resolution. 

Radio interferometers (such as ALMA) provide spa- 
tial and spectral resolution on interesting scales simul- 
taneously, making these instruments ideal observatories 
for probing substructure. We expect a significantly en- 
hanced sensitivity to substructure perturbations when we 
simultaneously model the visibilities observed in all the 
channels, compared to modeling the summed, channel- 



integrated visibility set. In the next section, we demon- 
strate this enhanced sensitivity using simulations of 
ALMA observations of lensed DSFGs. 



3. SIMULATIONS 

We simulate observations of lensed DSFGs with ALMA 
Cycle 1 (32 antennas) with the array in its most extended 
configuration, with a maximum baseline of 1.1 km. This 
allows a resolution of 0.16^' at 350 GHz. The noise levels 
are calculated using the online ALMA Sensitivity Calcu- 
lator for best observing conditions (1^^ Octile: PWV=0.5 
mm, Ts^s=90 K). Note that the fuh ALMA array wih 
have roughly 4 times greater sensitivity, and more than 
10 times higher spatial resolution than what we have as- 
sumed in the calculations presented here. We place the 
source at Dec -50:30 and simulate a one hour observa- 
tion with 10 second integration time. The uv-coverage 
was predicted using the "simdata" task of the Common 
Astronomy Software Application (CAS A) package (Petry 
et al. 2012). 



The source is modelled as a collection of giant star- 
forming clumps. Since we focus on high-excitation lines 
that are found only in high density clumps, we will ne- 
glect any diffuse emission originating from the interclump 
ISM. In local starburst galaxies, and in high redshift DS- 
FGs, CO line brightness peaks at J = 6^ and at higher 
J in AGN- dominated sources 
et al.|201Q[[Rangwala et al |20i 




as CO 6-5 produce some of the 
the dense cores of star forming clumps with compact 
morphologies, they may be the ideal molecular lines for 
substructure lensing detections. 

The clumps are placed randomly in both spatial loca- 
tion and velocity, with positions drawn from a gaussian 
profile with rms of 1 kpc and velocities drawn from a 
uniform distribution of width 300 km/s. Each clump 
has a circular Gaussian surface brightness profile with 
FWHM of dc = 250 pc, and a Gaussian profile in redshift 
space characterized by a velocity dispersion of cr^ = 30 
km/s. For simplicity we have used random velocities for 
the clumps, as opposed to a rotating disk, for example. 
As long as the clumps are reasonably well separated in 
velocity space (A'u ^ cr^) we do not expect any qualita- 
tive differences between ordered motion and disordered 
motion. Each source clump is therefore described by 6 
parameters: x and y centroids, central velocity, radius, 
flux and velocity dispersion. The source emission, sum- 
ming over all clumps, is represented as a 3D data cube. 
The total intrinsic velocity integrated fl ux of the entire 
DSFG is set to 1 Jy km/s (see BothwelTe TaT] (12012 ) 
for similar values for the velocity integrated CO flux of 
unlensed DSFGs) and distributed equally among 10 star 
forming clumps. 

Given our source realizations, we then ray-trace to 
compute a new 3D data cube containing the lensed im- 
ages observed in each channel. We then Fourier trans- 
form each layer of the cube, and predict the visibilities 
by interpolating the Fourier maps. 

In each simulation the main lens galaxy is modelled 
as a Singular Isothermal Ellipsoid (SIE), which has a 
three dimensional mass density proportional to r~^, and 
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Fig. 1. — Top Left: The source model and the lensed images of a clumpy source. The black curves show the tangential and radial 
caustics for the unperturbed macro model, and the red + symbol shows the location of an additional subhalo of mass M = 10^ M©. Top 
Right: The dirty image observed by ALMA. Bottom Left: Residuals of channel integrated dirty images between a smooth model and 
noisy perturbed observation. The greyscale is in units of noise rms. Bottom Right: Color residuals of the dirty images. The 50 observed 
channels are mapped to RGB colors as illustrated by the colorbar. The y-axis of the colorbar shows the intensity in units of image noise 
rms in each channel. 



projected surface density 



2G 



-1/2 



(1) 



where x and y are coordinates oriented along the princi- 
pal axes of the surface density measured relative to the 
lens centroid, v is the velocity dispersion along the line 
of sight, q is the axis ratio, and te is the Einstein radius 
of the lens. 



te = 47r- 



(2) 



( Kormann et al.|1994 ). The SIE therefore is described by 
5 parameters: the centroid Xc, ^c, velocity dispersion 



axis ratio and orientation Oq. 
Einstein radius mass Me = ttE 



We simulate lenses with 
cr / 1 of 4 X 10^^ Mq placed 
at Zd = 0.5 and place the source at z = 2. In addition 
to this main lens, we allow for external shear, described 
by an amplitude |7| and orientation 6^. The smooth 
mass model therefore has 7 parameters describing the 
lens. We follow previous work and model dark matter 
subhalos u sing the Pseudo-Jaffe density profile (Muhoz 
eraI][2QQT] ), 



1 

2x 



1 



(3) 



where x = r/rE and Xt = n/rE are the radius and the 



tidal truncation radius, respectively, in units of the sub- 
halo's Einstein radius, and = I](r)/I]cr is the dimen- 
sionless surface density expr essed in units of the lensing 
critical surface density (Sch neider et al.||1992 ). 

Figures [l] and [2] illustrate the increased sensitivity to 
substructure provided by spectroscopically resolved data. 
Figure [l] shows an example of ALMA observations of a 
lensed DSFG, both for spectroscopically resolved visibil- 
ities and channel-integrated visibilities. The greyscale 
and colored panels show the residuals from the best- 
fitting smooth lens models for each case, and the sub- 
structure perturbation clearly stands out more readily in 
velocity space. In Figure |2j we show the substructure 
parameter errors derived from simulated observations of 
another lensed DSFG, whose properties were chosen to 
be rep resentative of the lenses found in [Hezaveh et al. 
( 2Q12| ). The parameter uncertainties plotted in this Fig- 
ure were estimated by a Fisher matrix calculation; we 
marginalize over a considerably larger number of nui- 
sance parameters describing the source emission when 
fitting spectroscopically resolved visibilities, compared to 
fitting the channel-integrated visibilities. Despite the in- 
creased number of marginalized nuisance parameters, the 
parameter uncertainties are considerably reduced when 
we utilize the full, velocity-resolved data cube (compare 
grey vs. blue contours). This is true for the "macro 
model" parameters describing the smooth lens, and for 
the subhalo parameters as well. 
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Fig. 2. — Fisher matrix error forecasts. We simulate ALMA observations of a lensed DSFG, including substructure perturbations from a 
single subhalo of mass M = 10^ Mq in the lens mass distribution, and compute the Fisher matrix for all the lens (including substructure) 
and source parameters. We have marginalized over all parameters that are not shown in the figure. Gray contours correspond to fitting 
the channel-integrated flux, and blue contours correspond to fitting all channels simultaneously. 
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ALMA observations will therefore clearly improve our 
ability to characterize substructure in strong lenses. As 
Figure [2] demonstrates, the properties of a M = lO^M© 
subhalo can be determined with high precision, espe- 
cially when velocity information is utilized in the model 
fits. However, before undertaking a systematic study of 
subhalo parameter measurements, it is important to esti- 
mate the substructure detection rates expected in ALMA 
observations of typical lens systems. Will detections like 
the subhalo shown in Fig. [2] (a 6.Sa detection) be com- 
mon or rare in DSFG lenses? The focus of this paper 
is to demonstrate the benefits of velocity fitting, and to 
quantify the expected substructure detection rates using 
velocity fitting with ALMA. 

We define the significance of a subhalo detection in 
terms of the residual between the best-fit smooth 
mass model and the simulated observations. In our sim- 
ulations we start with a smooth macro model and add 
deflections from a single subhalo during the ray-trace. 
We then model the simulated observations using only 
the smooth mass model, with no substructure parame- 
ters, allowing the macro model parameters to adjust in 
order to acco unt for the substr uct ure perturbations. 

We follow |Dalal fc Kochanek| ( |2QQ2| ) and use linear 
perturbation theory to determine the parameters of the 
best-fitting smooth macro model. We characterize the 
smooth macro model with a parameter set p (including 
the source profile parameters). For some fiducial set of 
parameters po, we obtain observables O (a vector of real 
and imaginary visibilities) and write the penalty function 



X 



SO,{C-%SOj 



(4) 



where summation over repeated indices is understood, 
and C is the noise covariance matrix. If our observables 
are the visibilities, then C is diagonal. 

A perturbation to the model parameters Sp gener- 
ates perturbations to the observables given by SOi = 
{dOi/dpj)5pj. Given a current data residual SO^ the 
addition of these parameter adjustments changes the 
to 



X 



SOi + Spk 



dOj 
dpk 



(c- 



dpi 



Spi 



(5) 



To find the parameter set p that minimizes this func- 
tion we set dx^ /dpi = and solve for Spi 



Spi = -(F- 



dpk 



where 



F,,-— (C 



(6) 



(7) 



is the Fisher matrix for the smooth model. 

Using this procedure we choose a smooth macro model 
and calculate the matrix 80/ dp by finite differencing. 
This matrix has dimensions ?iobs ^ '^par? where Tiobs — 

18 

million is the number of observables (real and imaginary 
visibilities in 50 channels for a 1 hour observation), and 
^par = 57 is the number of source parameters (7 for 
the lens model, and 5 parameters eaclj^for 10 indepen- 

^ We did not vary each clump's velocity dispersion cr, since we 
did not expect it to be degenerate with any parameters of the lens 
mass model. 



dent source clumps). We use this matrix to calculate 
the Fisher matrix F using equation ^7 ) . We add a single 
subhalo to the current smooth model, predict the new 
visibilities, and use equation ^ to find the parameter 
adjustments to the smooth model that minimize x^ . We 
then simulate a smooth macro lens observation with pa- 
rameters p-\-5p and use the new x^ between this best-fit 
model and the original simulated observation as the like- 
lihood of a detection. 

We simulate four different lensing configurations la- 
beled "fold", "cusp", "double", and "cross". Ah configu- 
rations have parameters that are typical for strong lens- 
ing systems. As before, the background source is chosen 
to consist of 10 Gaussian clumps, of diameter (FWHM) 
250pc and velocity dispersion of 30 km/s. These values 
are chosen to be consistent with well-resolved observa- 
tions of molecula r gas clumps in lensed DSFGs ( [Swim] 
bank et al.||2010[ [2011) . Our sensitivity to substructure 
depends strongly on the sizes of these clumps, as dis- 
cussed in ^ For each configuration we define a Cartesian 
grid of substructure positions. We place a single subhalo 
at each location and use the above procedure to measure 
the detection significance of substructure with mass M 
as a function of position. The area inside a contour of 
specified x^ then defines the detection cross-section for 
each significance level. 

We perform simulations for 20 different substructure 
masses, between 10^ — lO^M©. Our results, shown in 
Figures [4] and |5j are discussed in the next section. 

4. RESULTS 

We define the substructure detection cross-section as 
the area of the sky (e.g. in square arcsec) inside of which 
a subhalo can be detected with a given minimum sig- 
nificance. We compute this cross-section by calculat- 
ing the detection significance, defined in ^ as a func- 
tion of subhalo position. Figure [3] shows the detection 
cross-sections at 3, 5, 7, and 10 sigma significance for a 
M = 10^ M0 subhalo, for four different macro lens con- 
figurations. The top row shows this cross-section when 
the channels are integrated before fitting and the bot- 
tom panel shows the cross-section when fitting to each of 
the channels individually. The different extent of the de- 
tection cross-section for different lensing configurations 
suggests that there is a higher probability of detecting 
sub-halos in high-magnification fold and cusp image con- 
figurations compared to low-magnification cross or dou- 
ble configurations. Nonetheless, the low-magnification 
lenses retain some sensitivity to substructure. This is 
unlike the case of 2-image quasar lenses, which generally 
lack sufficient constraints to permit substructure detec- 
tion. For DSFG lenses, the presence of multiple source 
clumps provides enough constraints on the mass model 
to permit substructure detection in favorable configura- 
tions. 

We perform similar simulations for other subhalo 
masses ranging between lO^M© to lO^M©, and calculate 
the total cross-section areas for each mass bin. Figure [4] 
shows these cross-sections as a function of mass for each 
lensing configuration. We have only plotted the cross- 
sections derived from channel- fitting, having already es- 
tablished that the detection sensitivity will be consider- 
ably diminished for channel-integrated data. 

These cross-sections, in combination with a substruc- 
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the model lensed images. The unlensed source is also plotted at the center. The black curves show the tangential and radial caustics, not 
the critical curves. The bottom row shows the detection significance when fitting to all channels simultaneously, while the top row shows 
detection significances for channel integrated visibilities. 




Fig. 4. — Detection cross-section for subhalos as a function of 
mass. Each panel corresponds to a different macro lens image con- 
figuration. The curves show the area of 3 (solid black), 5 (dashed 
blue), and 10 (dotted red) a detection significance. 
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Fig. 5. — Left: Cumulative mean number of detections of sub- 
halos of mass greater than M. The curves show the number of 
detections with 3 (solid black), 5 (dashed blue), and 10 (dotted 
red) sigma significance. 



ture mass function, can be used to estimate the ex- 
pected number of detected subhalos for each DSFG lens. 
We assume a subhalo mass function with a slope of 
(i log A/'/(i log M = —1.9, an upper mass bound of lO^M© 
and an overall normalization placed by setting the mass 
fraction of substructure inside the Einstein radius of the 
macro lens to / = 1%, consistent with results of the 



quasar lensing analysis of Dalai & Kochanek (2002). 
Multiplying the sky number density ot subhalos with 
masses between M and M + 5M by the detection cross- 
section of mass M gives the average number of detec- 
tions with the given significance, assuming a uniform dis- 
tribution of substructure over the strong lensing region. 
Figure [5] shows the cumulative number of detections for 
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the four lensing configurations. Note ttiat configurations 
witii tiigii magnification (tiie fold and cusp configura- 
tions) give more tiian 1 detected subiialo on average (at 
3(j confidence) for eacii lensing system. This somewhat 
invalidates our treatment, which explicitly assumes that 
only a single subhalo provides perturbations. In con- 
trast, the collective perturbations of multiple subhalos 
must be simultaneously treated. In forthcoming work, 
we will analyze the perturbations from the population of 
substructure inside dark matter halos. 

5. DISCUSSION 

In the simulations presented in this work, we have as- 
sumed a fixed source morphology with reasonably conser- 
vative parameters describing the DSFGs. In this section, 
we briefly explore the sensitivity of our results to the 
assumed properties of the source clumps. 

As discussed above, one of the most important pa- 
rameters describing the source is its angular size. In 
the calculations we have presented so far, we have as- 
sumed a physical size of 250 pc (FWHM) for the source 
clumps. We chose this valu e as a conservative upper 
bound, given the results of Swinbank et al. (2010). The 
CO-emitting regions of the star forming clumps could 
be significantly smaller than this assumed size, however. 
Upcoming ALMA observations of DSFGs should resolve 
the question of the intrinsic source size. Even if the 
source clumps are not significantly smaller than our as- 
sumed values, observations of other molecular lines be- 
sides CO lines could reveal mu ch smaller source sizes and 
different source morpholoffles ([Carilli et a L 2011 Riech- 
ers et al.|[2QTT| | Combes et al.||2 012'). I'ransitions ot H2U, 



HCN, and H(JO+ have higher critical densities than CO 
at similar observing frequencies. These lines trace dense 
gas in the active star forming regions (Downes fc Solomon 
[l998|), and may have fluxes as hi gh as 257o ot CO tlu x 



IE 



ao Solomon 2004) or more (Riechers et al. 2010). 
ecause they are contined to smaller regions than low and 
possibly high- J CO, they may increase the sensitivity 
to lower mass subhalos. To demonstrate this enhanced 
sensitivity, we have repeated the fold conflguration sim- 
ulations presented in Figure [4] with identical parameters, 
but using smaller source clumps (150 pc FWHM). The 
red dashed curve in the top panel of Figure [6] shows the 
5 a detection cross-sections of subhalos with masses rang- 
ing from 10^ to 10^ Mq. The blue shade corresponds to 
the larger clump sizes presented in Figure |4] The bottom 
panel shows the resulting increase in the mean number of 
detected subhalos, demonstrating that different molecu- 
lar lines originating from different morphological struc- 
tures of varying sizes exhibit different levels of sensitivity 
to substructure. 

Just as a smaller angular size improves the sensitiv- 
ity to substructure, narrower linewidths also help resolve 
the emission from distinct clumps. If the intrinsic ve- 
locity dispersion of each star forming clump is lowered, 
then fewer clumps appear simultaneously in each chan- 
nel. This effectively shrinks the angular extent of the 
emitting region in each channel, increasing the substruc- 
ture detection sensitivity. In Figure |6] the black solid 
curves show the 5<j detection cross-section (top) and 
mean cumulative detection numbers (bottom) of a DSFG 
composed of clumps with velocity dispersion of 10 km/s 
( [Hopkins et al.|2012 ). The key factor here is the relative 
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Fig. 6. — Detection sensitivity as a function of source properties. 
Here we show the improvement in detection significance when ob- 
serving DSFGs composed of smaller clumps (red dashed curves), 
clumps with lower velocity dispersion (solid black curves), or both 
(dot-dashed blue curves), compared to our fiducial source parame- 
ters (blue shaded region), for the fold configuration. The top panel 
shows the detection cross-section, and the bottom panel shows the 
mean cumulative number of detections. All curves are shown at 
5 a detection significance. Smaller clumps could correspond to 
observing dense-gas tracers such as HON. 



velocity offset of clumps with respect to each other, in 
units of their velocity dispersion. The blue dot-dashed 
curve in Figure |6] shows the detection improvements for 
simulations with both smaller clumps and lower velocity 
dispersion. As the bottom panel of Figure [6] illustrates, 
for a source with these parameters we expect to detect 
more than one subhalo in a given lens, if the substructure 
abundance is consistent with previous lensing analyses. 
Figure [6] demonstrates that our sensitivity to substruc- 
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ture depends on the clump morphology, both in real 
space and in velocity space. We have chosen fiducial 
values for these source parameters to be consistent with 
current measurements of the intrinsic proper ties of DS- 
FGs. The unlensed DSFGs presented in |Engel et al.| 
(2010) either showed signatures of merging systems, or 



fast rotating compact morphologies. Comparing the val- 



ues of galaxy - scale velocity dispersions reported iUjEngel 
et al I (|2010| ) which range between 200 and 1000 km/s 
(l^^WiilVl), wit h the ^ 10 km/s dis persions expected for 



single clumps (Hopkins et al. 2012) , indicates that all 
the DSFGs in their sample niust contain several com- 
ponents with significant velocity offsets. Similarly, the 
double peak velocity profile presented in [Hodge et al. 



(2012) is indicative of large velocity offsets between mul- 



tiple components. This suggests that systems with very 
small velocity offsets between their multiple components 
may be rare, so our assumption of Av > ay appears to be 
safe, even in the case of rotationally supported DSFGs. 

Our fiducial simulations used a fixed number of 10 
source clumps, randomly distributed across a galaxy of 
size 1 kpc. Our results do not appear to be as sensi- 
tive to these parameters as the clump size or linewidth, 
although this may depend on the details of the image 
morphology. For example, increasing the galaxy size al- 
lows the clumps to cover a larger portion of the caus- 
tic, which tightens the constraints on the macro model. 
However, it also increases the spacing between clumps, 
thereby degrading the sensitivity to low-mass substruc- 
ture. The effect of changing the number of clumps is 
also somewhat unclear. In principle, a larger number of 
clumps should provide a larger number of constraints. 
However, since we hold fixed the total flux, increasing 
the clump number makes each clump fainter and lowers 
their signal to noise ratio. We found that for the fold 
configuration described above, decreasing the number of 
clumps from 10 to 5 somewhat improves the detection 
sensitivity, mainly by making each clump brighter. It is 
unclear whether this will hold for the other configura- 
tions as well. In the absence of a systematic study all 
possible image configurations, we cannot say with any 
certainty how our substructure detection sensitivity will 
depend on the number and spread of the source clumps. 

Although in this work we have focused on high ex- 
citation molecular lines emitted by the compact cores 
of star forming clumps, the diffuse emission in low- 
excitation lines may also benefit from spatially resolved 
spectroscopy. Rotationally supported cold disks exhibit 
strong velocity gradients on the sky, meaning that the 
emission in narrow channels will originate from regions 
significantly smaller than the galaxy as a whole. In such 
cases the observed velocity gradients in the lensed im- 
ages of fast rotating cold gas reservoirs may show dips 
and peaks consistent with substructure lensing, comple- 
menting the high- J line observations. 

The calculations presented here were based on sim- 
ulations of 1-hour observations using ALMA Cycle 1. 
Longer observing times would improve our quoted sig- 



nal/noise ratios by t^/^. In addition, the simulations pre- 
sented here were specifically carried out for early science 
capabilities of ALMA with an array of only 32 antennas 
and a maximum baseline of ~1 km. The full ALMA array 
will consist of 64 antennas, which will increase the sensi- 
tivity by a factor of 4. In addition the longest baselines 
will extend as long as ~16 km, resulting in an angular 
resolution of ~10 milliarcsec at 850 /im. We can state 
with confidence that the full ALMA array will provide 
far greater substructure sensitivity than our Cycle 1 cal- 
culations have already found. 

In this work we have focused on ALMA observations 
of DSFGs, but in principle, the methods used here are 
not limited to interferometric observ ations. For exam- 
ple, integral field spec troscopy (e.g. Alaghband-Zadeh 



et al.|[2QT2{ [Barnabe et al. 2011; [ Forste r Schreiber et ah 
20111) of other lensed galaxies such Lyman Break Galax- 
ies (LBG) may benefit in the same way. 



6. CONCLUSION 

We have shown that it is possible to detect galactic 
dark matter substructure in the lens galaxies of lensed 
DSFG systems, using the early science capabilities of 
ALMA in Cycle 1. We found that the analysis of spatially 
resolved spectroscopic measurements of lensed sources 
can significantly improve substructure detection limits. 
In particular we simulated ALMA Cycle 1 observations 
of molecular lines in these systems, and calculated the 
detection significance for subhalos of different masses for 
various lensing configurations. We predict that with cur- 
rent ALMA capabilities and conservative assumptions 
about the morphology of DSFGs, there is a high proba- 
bility of detecting at least one subhalo with M > 10^ Mq 
in every lensed DSFG. Only marginally more optimistic 
source morphologies allow us to go below this limit and 
explore M ~ 10^ Mq dark matter subhalos. The full 
ALMA array, with 4x greater sensitivity and 16 x higher 
angular resolution than what we have assumed, will as- 
suredly improve our sensitivity to substructure beyond 
the calculations presented here. Given samples of ~100 
lensed DSFGs discovered in large area millimeter surveys 
and the significantly enhanced capabilities of the com- 
pleted ALMA, we are poised to revolutionize our under- 
standing of low-mass dark matter substructure in coming 
years. 
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